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Abstract 

We provide a complete calculation of all leading nucleon form factors at large momentum transfer 
and in single gauge boson exchange approximation. In order to evaluate the required nucleon 
transition probability matrix elements, we combine QCD perturbation theory with an expansion 
in nucleon distribution amplitudes. Using leading twist nucleon distribution amplitudes only, one 
obtains the desired nucleon form factors. The obtained results are consistent with experimental data 
when we use third order polynomials for the distribution amplitudes and non-perturbative models 
for the coefficients. Finally, we compare our results for the Dirac form factor of the proton with 
previous results based on the QCD factorization theorem, where a discrepancy about the correct 
symmetry factor is still under discussion. We examine the calculations and propose an explanation. 
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I. INTRODUCTION 



Studies about leading nucleon form factors received much attention in the past. In the 
region of large momentum transfer, the formulation of the QCD factorization theorem in 
|lL fl, fl, and Q|, js| was a great success. In [2] and [3|, the expressions for the leading 
electromagnetic form factors were derived. Using the QCD factorization theorem, these 



expressions were improved in 



7]. Moreover, the representation for the isovector axial- vector 



form factor was created in [8] with applications discussed in Similarly, the isoscalar 
axial- vector form factor was considered in [10]. Studying the behavior of vacuum to nucleon 
projection matrix elements, the leading electromagnetic form factors were calculated in {ll| 



and [12] at the asym 
was also applied in 



atotic limit. Hereby, an averaged value of a s was used. This approach 
13 1 and improved in fl^ij . Furthermore, in {l4|, also non-asymptotic 



contributions of the leading twist nucleon distribution amplitude were taken into account 



to calculate the leading electromagnetic 
keep a s inside the integral, evaluated in 



formalism was criticized in 



brm factors, see also 



15]. Another approach is to 



161 ] . Despite the success of perturbative QCD, the 



17fl. Therefore, a modified formalism was discussed in [18|. Being 



more precise, the region of low momentum transfer was considered in the case of the pion, 



see 
see 



19 1 . Using the modified formalism in the case of the pion, the problems were cleared, 
20|. This technique was also applied for the nucleon in j2l|. An improved calculation of 



the proton magnetic form factor, within the modified factorization scheme, was given in 22]. 
An analogous calculation for the neutron was done in 23J. A comprehensive work about 

n n 

the theory was published in [2J|. A remaining problem has been pointed out in [25], where 



a different symmetry factor was used as in 



16|. The same difference was also discussed in 



26|. Meanwhile, nucleon distribution amplitudes of sub- leading twist were also extensively 



studied, see |27fl, with applications in [28]. 

This work is organized as follows. First, we will introduce the general construction 
for the desired form factors. Second, we will calculate the leading nucleon form factors. 
Next, we will compare our results with experimental data. Afterwards, we will discuss the 
normalization factor discrepancy. Finally, we draw our conclusions. Important technical 
details are collected in three appendices. 
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II. GENERAL CONSTRUCTION OF THE NUCLEON FORM FACTORS 



Elastic lepton nucleoli scattering in the single gauge boson exchange approximation is the 
basic process to receive information about the nucleon structure within QCD. Calculating 
the cross section of these processes, one has to study nucleon transition probability matrix 
elements including a single quark current. It is well known that they can be expressed 
by nucleon form factors. Expanding these matrix elements at large momentum transfer 
Q 2 = —q 2 = — (P' — P) 2 , where P and P' are, respectively, the incoming and outgoing 
nucleon momenta, one can reduce the decomposition to one dominant term depending on 
the leading nucleon form factor. Overall, one gets four independent expressions. 

It is noteworthy that the computation of the different form factors is similar from the 
theoretical perspective, whereas they differ in the experimental realization. 

Nucleon matrix elements including the electromagnetic current can be expressed in terms 
of the the magnetic form factors or, equivalently, the Dirac form factors in the specified case, 
using obvious abbreviations, 

(p(P0|JT(O)b(P)> = G p M (Q 2 )N(P'h,N(P) (1) 
(n(P')\J-(0)\n(P)) = G n M (Q 2 )N(P'h,N(P). (2) 

These structures describe the exchange of a highly virtual photon that is the dominant 
process in experiments. We notice that they are also part of neutral weak interactions. 

In the the axial sector, we use the expansion in terms of the axial- vector form factors. In 
order to calculate them, it is necessary to choose the processes 

(iV/(P / )l<(0)|iV l (P)) = lG^Q 2 )N(P>) w 5 fi N(P) (3) 

(N f (P')\A«(0m(P)) = \g v a (Q 2 )N(P') W t^N(P)- (4) 
We work with two proton states for G S A , meaning that we have to evaluate the matrix el- 
ement (p(P') 1^(0)7^5^(0) +^(0) 7m75 ^(0)|p(P)> = G A (Q 2 )N(P') llMl5 N(P). Concern- 
ing G V A , we consider the proton to neutron transition. We end up with the reduced matrix 
element (n(P') 1^(0)7^(0) b(P)> = G A (Q 2 )N(P')^ l5 N(P). 

The isovector current matrix elements describe a part of the charged weak interactions 
generated by a W + or W~ boson and the neutral weak interactions mediated by a Z° boson. 
The isoscalar current matrix elements only appear in non standard weak interaction theories, 
such as the exchange of an extra Z° boson. 
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III. CALCULATION OF THE LEADING NUCLEON FORM FACTORS 



Starting from the given four matrix elements, we make use of the large momentum transfer 
for a perturbative expansion. Therefore, we need the S-Matrix containing the interaction 
part of the QCD Lagrangian. In all cases, one has to expand up to the fourth order to 
obtain non-vanishing results. We mentioned that the calculations of the different form 
factors are similar, so we demonstrate the evaluation of the most prominent member, the 
proton magnetic form factor 



\P(P))- (5) 



This expansion can be expressed by Wick Contractions and Feynman Diagrams. These 
diagrams have similar structures and therefore, we will calculate just one diagram in detail. 
The entire diagram classification will be specified in the appendices. 



A. Sample Diagram Evaluation 

The chosen diagram can be presented as follows. We have the incoming proton on the 
left and we have the outgoing proton on the right. The cross denotes the coupling to the 
virtual photon and the quark lines denote u, u, d from top to bottom. The designations 
at the vertices are the corresponding coordinates and the designations at the lines are the 
corresponding momenta. 




— »> * 

u 3 P x * v 3 P' 



Let us begin with the determination of the color factor. Therefore, one has to examine 
the color structure of the diagram, denoting the color indices with (a, . . . 
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i i i i i i i i 

(p(POI[^l)]&[4(z 3 )]/[^ 

Combining all terms and contracting the generators, one gets the result of the color factor 

This result is the same for all diagrams of G P M , G 7 ^ and G S A . The flavor change in the 
diagrams of G V A leads to the exchange of e^fh — > £bhf producing an overall minus sign. 

We continue with the evaluation of the Lorentz structure of the diagram, designating 
the Lorentz indices by (a, . . . One gets an expression as part of (jSJ), where we already 
included so that 



''^■ e «n / ^ 4x ihp]ab[l ai }cd[la 2 }ef[l a3 }ghh 



54 

i=l 



[Mxi)UM0)}a [^ u {x 3 )} h [^ u {x 2 )} e A^{x 1 )A^{x 2 )A^{x 3 )A ai {x i ) 

I I I I I - I I I I 

(p(POI$«(*l)]c0«(z 3 )]^ 

The formulas for the propagators are standard and the representations of the vacuum 
to nucleon projection matrix elements are extensively studied in [27]. Applications are 
discussed in [28|. Hereby, g aiOL i is the metric tensor. Recombining the Lorentz structures, 
one gets the following result 

i=i j v ' j=i j j k=\ j k j 

g— ia;i-(Ai-Ai— v\p') ^-ix2-{-^2+^i+U2p) ^-ixz-{/^2-^2-V2p') ^-ixi-(k2+uzp-vzp') 

The integration over the coordinates produce delta functions which can be used to inte- 
grate over the momenta. They describe the momentum conservation at each vertex of the 
diagram 

Aj. = (v i + v 2 + v 3 )p' - (u 2 + u 3 )p A 2 = (« 2 + v 3 )p' - u 3 p 

Ai = (v 2 + v 3 )p' - (u 2 + u 3 )p A 2 = v 3 p' - u 3 p. 

The undeclared component S is the sum of all remaining structures connected with 
combinations of nucleon distribution amplitudes, where V, A, and T denote vector, axial- 
vector, and tensor structures, respectively, and nucleon spinors. For convenience, we omit 



the twist index at the distribution amplitudes and the dependence on the quark momentum 
fractions. Moreover, we only have to deal with the large component of the spinor, so we can 
use the standard notation. One gets 

S, = N(P')y a tN(P)Tr[y l jfy a Jt 2 y a Jy Ql £ 1 )(VV + AA) 

5 2 = N(P')y a4 y 5 N(P)Tr[y^y a ^ 2 y a3 y^y ai ^ 1 ](-AV-VA) 

5 3 = iV(P / )7 A '7a 4 7 A iV(P)Tr[7^cT Ap 7a 2 A27 a 3^v P '7a 1 Ai](-rT). 

Finally, we can express the result of the discussed diagram depending on the integration 
over the quark momentum fractions in the form 



(4:7ia s 



[du] 



[dv] 



216 Q 4 7 u 3 (u 2 + u 3 ) 2 v 3 (v 2 + v 3 f 



[(V - A) 2 + AT 2 } N(P')y^N(P). 



(7) 



B. Results for the Form Factors 



We add the obtained results of the contributing diagrams for every form factor and insert 
the pol ynomial expression of the independent leading twist nucleon distribution amplitude, 
see 24] and J27 1. As already explained, we apply third order polynomials expressible as 
$3(^1, u 2 , u 3 ) = 120uiU 2 u 3 fN[ci + c 2 Ui + c 3 u 3 + c^ul + C5UI + CQU1U3}. The required coefficients 
are non-perturbative parameters which are scale dependent and have to be modeled or 
taken from the lattice. We use the standard hadronic scale of 1 GeV, and adopt /jy = 
5.3 ±0.5) ■ 10~ 3 GeV 2 , see 27j. Meanwhile, this parameter is also studied on the lattice, see 
29 1 . The first and second order coefficients are taken from {27 1 and the third order modeled 



coefficients can be found in [24]. We use the models of Chernyak, Zhitnitsky (CZ) [111 ]. Gari, 
Stefanis (GS) (l3|, King, Sachraida (KS) H, Chernyak, Ogloblin, Zhitnitsky (COZ) jl5|, 
and Stefanis, Bergmann (HET) |3jJ. The abbreviation HET means heterotic conception. 

Every form factor can be expressed by a corresponding coefficient function as Q i G{Q 2 ) = 
(4ira s ) 2 ffjX(ci, c 2 , c 3 , C4, C5, Cq). Finally, we have the following list of them: 



X p M (c x , c 2 , c 3 , c 4 , c 5 , c 6 ) = Yl Yl a v CiC i ( 8 ) 

i=l j=l 

with a u = 0, a 21 = 4320, a 3l = -1728, a A1 = 6048, a 5 i = 432, a 61 = -1296, a 22 = 2 1 60, 
a 32 = 672, a 42 = 4536, a 52 = 1152, a 62 = -240, a 33 = -528, a 43 = 1464, a 53 = -168, 
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= —480, Gt44 = 2123, 054 = 1418, = 72, a 55 = 95, = —180, = —72, 



Xm(ci, c 2 , c 3 , c 4 , c 5 , c 6 ) = ^2 a ij c i°j ( 9 ) 

i=i j=i 

with an = 3888, 021 = 432, 031 = 3456, 041 = —972, 051 = 1836, a 6i = 648, 022 = —336, 
a 32 = 144, a 42 = -960, a 52 = -180, a 62 = 120, a 33 = 1008, a 43 = -336, a 53 = 1392, 
a 63 = 240, a 44 = —491, 054 = —395, a 64 = 0, 055 = 523, = 126, a 66 = 24, 



-^l( C l5 C 2; c 3 5 C45 C 5 , C 6 ) — — — ^ ^ dijCiCj (10) 

i=i i=i 

with an = 1296, a 2 i = 1872, a 3 i = 144, a 4i = 1872, a 51 = 216, a 61 = -144, a 22 = 640, 
a 32 = 448, a 42 = 1220, a 52 = 432, a 62 = -8, a 33 = -176, a 43 = 512, a 53 = -25 2, a 63 = -72, 
0,44 = 549, a 54 = 440, a 64 = 42, 055 = —100, = —30, = —12, 



-^a( c i> c 2; c 3 ? C4, c 5 , c 6 ) — — gggg ^2 ai J CiC i (1-0 

1=1 jr = l 

with an = 7776, a 2 i = 3456, a 3 i = 3456, 041 = —216, a.51 = —216, a 6 i = 1728, 022 = —304, 
a 32 = 1472, a 42 = -1560, a 52 = 312, a 62 = 480, a 33 = -304, a 43 = 312, a 53 = -1560, 
a 63 = 480, a 44 = -907, a 54 = -124, a 64 = 144, a 55 = -907, a 65 = 144, a 66 = 72. 



IV. COMPARISON WITH EXPERIMENTAL DATA 



The final step is to compare the results for the different form factors with available 
experimental data. In order to compare our results with those existing in the literature, we 
do not specify the values of the parameters /jy and a s . 

The determination of a s has to be done with care, because it depends on the momentum 
transfer of the process, or more precisely, it depends on the gluon virtualities. Consequently, 
one gets two of them for every diagram independently. The main problem arises from the 
soft gluon region, where a s becomes divergent. It is possible to apply different estimations 
;o avoid this problem. In principle, one can introduce an averaged a , fo r all diagrams, see 
11 1, jl2 1 . This approach was also applied in 13 1 and improved in 14 1 for each diagram 
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separately. Hereby, the value of a s ~ 0.3 was used. Another approach is to keep a s inside 



the integral, as done in 



16] . This technique can change the slope of the form factors and 



improve the scaling behavior of them. Nevertheless, this cannot be done without introducing 



a new cutoff parameter. Taking into account soft gluon corrections [21| . one may claim that 
the area of soft gluon virtualities is systematically underestimated. According to this, we 
prefer to use an averaged value of a s which is taken from an effective momentum transfer 
around 1 GeV 2 . Consequently, we use a s = 0.45 ± 0.05. One should note that the obtained 
results contain large uncertainties due to the choice of a s . 

The first order polynomial describes the asymptotic limit of Q 2 leading to the following 
results Q A G P M {Q 2 ) = 0(4vra s ) 2 /^ = OGeV 4 , and Q 4 G n M (Q 2 ) = 50/3(4vra s ) 2 / 2 = (15 ± 
5) ■ 10- 3 GeV 4 , and Q A G V A {Q 2 ) = 50/3(47ra s ) 2 / 2 r = (15 ± 5) • lO" 3 GeV 4 , and Q A G S A {Q 2 ) = 
-50{4na s ) 2 f 2 = -(45 ± 5) • 10" 3 GeV 4 . 

The second order polynomial is still related to values of Q 2 which are unreachable in 
experiments, but one obtains a non-asymptotic behavior. We get the results Q 4 G^/(Q 2 ) = 
lQQ(Aira s ) 2 f^ = (0.15±0.05) GeV 4 , and Q 4 G^(Q 2 ) = -74(4vra s ) 2 /^ = -(0.07±0.05) GeV 4 , 
and Q 4 G V A {Q 2 ) = 247(4vra s ) 2 / 2 r = (0.22 ± 0.05) GeV 4 , and Q 4 G S A {Q 2 ) = 256(4vra s ) 2 / 2 r = 
(0.23 ± 0.05) GeV 4 . 

The third order polynomial can be correlated to values of Q 2 in the area around 10 GeV 2 . 
In this case, one can compare with experimental data, though one should keep in mind that 
the soft contribution may still be sizeable, see 22j and 



For the proton magnetic form factor we obtain 



Q 4 G P M {Q 2 ) = (47ra,) 2 / 2 ■ < 



1.16 


10 3 




' (1.0 ±0.1) 


GeV 4 


(CZ) 


1.16 


10 3 




(1.0 ±0.1) 


GeV 4 


(GS) 


1.69 


10 3 


> = < 


(1.5 ±0.1) 


GeV 4 


(KS) 


1.34 


10 3 




(1.2 ±0.1) 


GeV 4 


(COZ) 


1.53 


10 3 




k (1.4 ±0.1) 


GeV 4 


(HET) 



(12) 



The electromagnetic interaction ensures that this form factor can be measured at large Q 2 . 
Moreover, the proton is a stable particle and so the measurement can directly be done on 
the proton. We have data of G P M for a comprehensive region, see [32] and j33|. The obtained 
results are in sufficient agreement with the averaged experimental value Q 4 G P M (Q 2 ) = (1.0 ± 
0.1) GeV 4 . 
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For the neutron magnetic form factor we find 

' -0.56 • 10 3 
-0.11 • 10 3 

Q 4 G n M (Q 2 ) = {^a s ff 2 N 



-(0.5 ±0.1) GeV 4 (CZ) 
-(0.1 ±0.1) GeV 4 (GS) 
-(0.6 ±0.1) GeV 4 (KS) 
-(0.6 ±0.1) GeV 4 (COZ) 
-(0.1 ±0.1) GeV 4 (HET). 



(13) 



-0.70 • 10 3 
-0.63 • 10 3 
-0.16 • 10 3 

The electromagnetic interaction ensures that this form factor can also be measured at large 
Q 2 . Unfortunately, the neutron is not a stable particle and so the measurement cannot 
be done on the neutron directly. Data of G^j in the lower region of Q 2 , see [34] and 
|35j, are available, but an extrapolation to the upper area is also possible. The obtained 
results are in sufficient agreement with the extrapolated and averaged experimental value 
Q 4 G^(Q 2 ) = -(0.5±0.1)GeV 4 . 

The next result concerns the isovector axial- vector form factor 

1.6 ±0.1) GeV 4 



Q G V A {Q ) = (4vra s )^ 



' 1.75- 10 3 ' 




1.31 • 10 3 




< 2.50 • 10 3 


> = < 


2.04- 10 3 




1.80 • 10 3 

V J 





(1.2 ±0.1) GeV 4 

(2.2 ±0.1) GeV 4 

(1.8 ±0.1) GeV 4 (COZ) 

H a _i_ n 1 \ „t r4 /ttt^t^ 



(CZ) 
(GS) 
(KS) 



(14) 



This form factor appears in standard weak interactions. The consequence is that the desired 
processes are generally suppressed by the dominant electromagnetic interaction, but it is 
possible to choose processes for the measurement which can be separated from the electro- 
magnetic interaction. Avoiding the electromagnetic domination, one has to apply neutrino 
neutron scattering. The neutrino interacts slightly and so we have data for values of Q 2 
which are less than required, see 36]. Predicting the behavior in the required region, one 
can extract a value for G\ which is given by Q 4 G V A (Q 2 ) = (1.5 ± 0.2) GeV . 
Last, we consider the isoscalar axial- vector form factor 



Q*G A {Q 2 ) = (4na s ) 2 f 2 N 



' 1.66 


10 3 ' 




' (1.5 ±0.1) GeV 4 


(CZ) 


3.10 


10 3 




(2.8 ±0.1) GeV 4 


(GS) 


< 2.61 


10 3 


> = < 


(2.3 ±0.1) GeV 4 


(KS) 


1.89 


10 3 




(1.7 ±0.1) GeV 4 


(COZ) 


3.71 

V 


10 3 

J 




k (3.3 ±0.1) GeV 4 


(HET) 



(15) 



This form factor appears in non standard weak interaction theories only. Consequently, the 
required processes are suppressed by all discussed interaction theories. Therefore, it is not 
possible to construct an experiment in order to measure this form factor independently. We 
have experimental data on G S A for values of Q 2 which are significantly less than required. 
Under these circumstances, it makes no sense to compare the results. According to this, we 
can only consider our results as a prediction for further experiments. 



V. DISCUSSION OF THE NORMALIZATION FACTOR DISCREPANCY 



Finally, we will compare our results for the proton magnetic form factor with previous 



results obtained in calculations based on the QCD factorization theorem, see 



1] 



and 



flh Q|j S ■ Hereby, we focus on the known discrepancy concerning the symmetry factor. The 
starting point of our considerations is the calculation of the desired form factor performed 
in {if]]. After some time, the discussed form factor was recalculated in {25 1. They obtained a 
result which is precisely a factor of 2 smaller than the previous result. They mentioned that 
they do not understand the origin of this discrepancy. A newer calculation of the proton 
magnetic form factor was done in (26 1. Their result is in agreement with (25I . They realized 
that the symmetry factor in is precisely a factor of y/2 larger than in (25] - Nevertheless, 
they mentioned that they do not know which symmetry factor is correct. However, our 
result of the proton magnetic form factor is the same as in 26j. At this point, we should tr 



noma try 
ing in (3|, 



to find the origin of the discrepancy between the older and the newer results. Looking 
we can see that the normalization of the soft parts is not fixed, but a required normalization 
connected to the number of colors is discussed. We realized that we need a normalization 
connected to the number of possible color states. In the case of the pion, one gets no 
difference, because the normalization by the number of colors as well as the normalization 
by the number of possible color states would lead to the factor 1 / y/3. This is not true in 
the case of the nucleon. The normalization by the number of colors would obviously lead 
to the factor l/\/3 too, but the normalization by the number of possible color states would 
require the factor l/v6 creating a difference of y/2. Because two soft parts are required, 
one has to square the normalization factor. We claim that this is the source for the factor 



of 2 difference in previous calculations. Indeed, when we compare the normalization of 
and [25] , we find the normalization factor 1 / y/3 in 16| and 1 / v6 in [25] as expected. 



16| 
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VI. CONCLUSION 



In this work, we calculated all leading nucleon form factors at large momentum transfer 
and in single gauge boson exchange approximation. 

In order to compare the obtained results with experimental data, we used third order 
polynomials for the distribution amplitudes and non-perturbative models for the coefficients. 

The calculated results for the form factors are essentially in good agreement with the 
available experimental data. Nevertheless, further experimental data at large momentum 
transfer are required in order to gain a deeper understanding of the nucleon form factors. 

Finally, we discussed the known discrepancy concerning the symmetry factor in previ- 
ous works. We argued that taking into account the color normalization appropriately, the 
discrepancy can be resolved. 



Appendix A: Diagrams 

The general decomposition in diagrams can be classified by the quark coupling to the 
specified current in the corresponding matrix element. In principle, one gets 14 diagrams 
for every quark. The only exception is that the transition expression of G\ does not allow 
the coupling to the quark which occurs once in the initial nucleon. The discussed diagrams 
are not independent, because two quarks in the initial nucleon are identical and a general 
interchange of the incoming and outgoing nucleon does not change the final result. 



Fixing the positions of the quark lines in the manner that the single quark in the initial 
nucleon is at the top or bottom position, one has to deal with 7 diagrams as follows. 




We can derive the following list of momentum conservation constraints when the single 
quark in the initial nucleon is at the bottom position. The expressions for the diagrams with 
the considered quark at the top position can be obtained by the interchange of u\ -H- u 3 
and v 1 <->■ v 3 . This list can be used for all form factors except for G V A where an exchange of 
v 2 <-> v 3 is required 

Ai = (vi + v 2 + v 3 )p' - (u 2 + u 3 )p A 2 = (v 2 + v 3 )p' - u 3 p 

Ai = (v 2 + v 3 )p' - (u 2 + u 3 )p A 2 = v 3 p' - u 3 p 
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1 — 1 ± 


= (v-i + 1>9 + v-i)v' — (llo + u^)v 

\ J- 1 £ 1 O/Jr \ Z 1 o 


A 2 


= (ur> + Ui)v 

\ ^Z i o J ir 


— V-xX)' 


Ai 




A 2 


= v 3 p' - u 3 p 




1 — 1 ± 


= (v-i -\- Vr> -\- V-i)t>' — (ur> + U^)v 


A 2 

1 — "Z 


= (v<) + v^)v' 




Ai 


= (v 2 + v 3 )p' - (u 2 + u 3 )p 


A 2 


= u 2 p - v 2 p' 




Ai 

1 — * 1 


= (vi + + V'?)??' — (it') + wO/D 

V 1 1 ^ Z 1 ^ o } \^Z 1 ^O/r 


A 2 

1 — *z 


= (it') + "U^/D 


— Vol/ 


Ai 


= (V 2 + V 3 )p' - (ti 2 + «3)P 


A 2 


= W 2 £> - V 2 p' 




Ai 

1 — * 1 


= (v~t -\- Vo -\- V-i)!)' — (ur) + U-\)v 
\ 1- 1 ^ 1 o J Jr \ z i ^o/lr 


A 2 

1 — *z 


= ll>i + VilT)' 




Ai 


= v 2 p' - u 2 p 


A 2 


= v 3 p' - u 3 p 




Ai 


= (vi +v 2 + v 3 )p' - (u 2 + u 3 )p 


A 2 


= (vi + v 2 )p' 


- u 2 p 


Ai 


= v 3 p' - u 3 p 


A 2 


= v 2 p' - u 2 p 




Ai 


= (ui + u 2 )p- v 2 p' 


A 2 


= (vi + v 3 )p' 


-u 3 p 


Ai 


= v 2 p' - u 2 p 


A 2 


= v 3 p' - u 3 p. 





Appendix B: Structures 

Here we present the structures summarized in the introduced element S. These compo- 
nents are similar for all form factors except for G V A . 

When the single quark in the initial nucleon is at the bottom position, we get the following 
list of structures for all remaining form factors except for G S A . Its results require the exchange 
of 7^ -» 7 M 7 5 . We find 

5 1 = N(P')^ a Mn^h^a 2 ^aJ^i}(VV + AA) 

5 2 = 7V(P')7a 4 75^(P)Tr[7>7 a2 A27a 3 75/7a 1 Ai](-^-^) 

5 3 = iV(P')7 A 7o4 7 

N(P) Tr[7 M io- Ap 7 a2 ^ 2 7 a3 ior A / p /7 ai ^: 1 ] (-TT) 

5 1 = N{P') lai N{P) Tr [ 7 >7a 3 A27a 3 /7a 1 ^(VV + AA) 

5 2 = ^(P07a 4 75^(^)Tr[ 7 >7 a3 7 5 A27ay7a 1 Ai](-^-^) 

5 3 = iV(P / )7 A '7a47 A A r ( J P)Tr[7^ ( TA P 7a3A27 a2 ^v P '7a 1 Ai](-TT) 
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Si 

s 2 

S3 

Si 

s 2 

S3 

s 2 

S3 

s l 
s 2 

S3 

S! 

s 2 

S3 



N(P' 
N(P' 
N(P' 

N(P' 
N(P' 
N(P' 

N(P' 
N(P' 
N(P' 

N(P' 
N(P' 
N(P' 

N(P' 
N(P' 
N(P' 



lai /K 2la2 N{P) ^\ lv tf la J lax PL x \{VV + AA) 
l ai ls/K 2la2 N(P) Trfr^^Mi-AV - V A) 

I012I 

N(P) Tr[y f ,i(T Xp 'y a:i i(Txp>'y ai jXi\ (-TT) 

la2 /K 2lai N{P) Tr [7^3/7^1]^ + AA) 
la2 /K 2lailb N(P) Tr [i^l^la, £i](-AV - V A) 

7a 4 7 

lai N{P) Tr[ 7A y 7a y 7^27^1]^ + AA) 
lailb N(P) Tr[7 A y7 a y7 a3 7 5 A27a 1 Ai](-^^ - V A) 
7 A 'la A l X N(P) Tr[7 M ia Ap 7 7a3 A 27Q1 Ai](-TT) 

7a2 iv(p) Tv[7>7 a y 7Q3 A 2 7 ai A 1 ](yy + aa) 

7a2 75iV(P) Tr[ 7A y7 a y7 a3 A27a 1 75Ai](-^^ - VA) 
7 A '7a 2 7 A ^(-P) Tr[7 M ia- Ap 7 

7asA27aiAl] (-TT) 

lai N{P) Tr[7 M Ai7a 1 /7 a y7a 3 A2](^V + AA) 
7a4 7 5 iV(P) Tr[7 M Ai7ay7ay7a 3 75A 2 ](-A\/ - IM) 
7 A ' 7a 4 7 A ^ (P) Tr [7^1 7ai ^Ap7a 2 ItTAV 7a 3 A2] ( -TT) . 



When the considered quark is at the top position, we get again a list of structures for all 
remaining form factors except for G S A . Its results require the same exchange as applied in 
the previous case. We get 

5! = iV(P / )7 ai Ai7M^(^)Trb7 a2 A27 a y7 a4 ](^V + AA) 

5 2 = iV(P / )7 ai 75Ai7M^(P)Tr[^7 a2 75A27 a y7 a4 ](-A\/-\/A) 

5 3 = iV(P / )7 A '7a 1 Ai7M7 A iV(P)Tr[m Ap 7 a2 A27a3^v P '7 a4 ](-TT) 



5 1 = N(P') laiJ X 1 ^N(P)Tr\^ a3 jX 2la J la A(VV + AA) 

5 2 = 7V(P')7a 1 75Ai7^(P)Tr^7 Q3 A 2 7a 2 7y7a 4 ](-Ay-yA) 
5s = iV(P')7 A '7 ai Ai7M7 A iV(P)Tr[ ? a Ap 7 a 3A27 

i(y\' P 'la 4 ]{-TT) 
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51 = N(P')^ ai ^N(P)Tr\^ a ^ a4 ^ a2 ](VV + AA) 

5 2 = NiP^w^^N^TrW^^wM-AV-VA) 

5 3 = iV(P / )7 A '7a 1 Ai7M7 A iV( J P)Tr[ ? a Ap 7 a 3 ?C T AV 7 a4 A27 a2 ](-TT) 

51 = iV(P / )7 ai Ai7^( J P)Tr^7 a y7a 2 A27 a4 ](^^ + AA) 

5 2 = iy(P')7 ai 75Ai7^(P)Tr^7 a y7a 2 75A27 a 4](-^^-\/A) 

5 3 = iV(P')7 A '7a 1 Ai7M7 A iV(P)Tr[ ?( T Ap 7 Q 32 ( T AV 7 a2 A27 a4 ](-TT) 

51 = N(P')^ a3 ^ ai ^N(P)Tr\^ a ^ a4 ](VV + AA) 

5 2 = 7V(P07 Q3 A27a 1 75Ai7^(P)IV^7 a2 7 5 ^7 Q4 ](-Ay-yA) 

5 3 = iV(P / )7 A '7a3A27a 1 Al7/,7 A ^(P)Tr[2(T Ap 7a 2 ^V P '7a4](-^) 

5 X = N(P , )^ 2 ^ ai ^N(P)Tr]^^ m ](VV + AA) 

5 2 = iV(P07a 3 A27a 1 75Ai7^(P)Tr^7 Q y7 a2 7 5 ](-A\/-\/A) 

5 3 = iy(P / )7 A '7a 3 A27a 1 Ai7M7 A ^(P)Tr^ Ap 7 Q4 2av p ,7a 2 ](-TT) 

51 = iV(P')7 a3 A27, t Ai7a 1 iV(P)Tr^7 Q y7 a4 ](^^ + AA) 

5 2 = iy(P')7 a3 A27, t Ai7 Ql 75iV(P)Tr[/7a 2 75/7 a4 ](-^^-\/A) 

5 3 = iV(P / )7 A '7a 3 A27^Ai7a 1 7 A ^(P)Tr[2(T Ap 7a 2 ^v p ,7a 4 ](-TT). 

We obtain the following list of structures for G\. The discussed quark must be at the 
bottom position here, so that 

51 = N(P') la3 /X 2la ^ l5 ^ aa ^ la4 N(P)(VV + AA + AV + VA) 

5 2 = iV(P / )7 A '7a 3 A27ay7M75Ai7 ai ^Av7a 4 iV(P)(-T\/ + TA) 

5 3 = iV(P')7 a3 A27a 2 ^A P 7M75Ai7ay7 a4 7 A iV(P)(VT-AT) 

5 4 = iV(P , )7 A '7a 3 A27a 2 ^A P 7M75Ai7a 1 ^v P '7 a4 7 A ^(P)(-TT) 

51 = N(P , )l a JC2l a ^tf5£ll«jflo«N(P)(VV + AA + AV + VA) 

5 2 = iV(P / )7 A '7a 2 A27a 3 y7M75Ai7 ai ^v P '7a 4 A r (^)(-TV + TA) 

5 3 = iV(P')7 a2 A27a 3 ^A P 7M75Ai7ay7 a4 7 A iV(P)(yT-AT) 

5 4 = iV(P / )7 A '7a 2 A27a 3 ^A P 7M75Ai7a 1 ^A' P '7 a4 7 A iV( J P)(-TT) 
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Si 




la3 ^^ lla J laA /K 2la2 N{P){VV + AA + AV + VA) 


s 2 


= iV(P'' 


A' / A/ * \y at/ 7 ~> \ / mT t 1 m A \ 
7 7a.-y7^75Al7a 1 ^A'p'7a 4 A27a 2 A^(P)(-T\/ + 274) 


S3 




7 a >A»7u7 5 Ai 7^/7^27^7^^) (^T - AT) 


s 4 


= 


7 A '7 Q ,JaAp7»75Ai7ai^ "A'p'7a 4 A27a 2 7 A iV(P)(-PP) 


<Si 


= 


l a jfl^/K lla ^ lai /K 2la2 N{P){VV + AA + AV + VA) 




= N(P r 


A ' X A/ * A/ A 7" / 7~ > \ / m T 7" 1 m A \ 

7 7ay7^75Ai7a3^A' P '7 a4 A27a 2 ^(P)(-T\/ + 774) 


s 3 


= N(P'] 


7 ai ^A P 7u7 5 Ai7 a y7a4A27 a2 7 A A r (^)(^ / T - AT) 




= N(P'] 


7 A '7ai^ "Ap7/ t 75Al7a3^ "A'p'7a 4 A27a 2 7 A iV(P)(-TT) 




= N{P'] 


la jf llllb /K llai /K 2la J lai N{P){VV + AA + AV + VA) 


S 2 


= N(P r 


A' / A/ AX * AT/ T~l\ / /"7~T T T 1 /Tl \ 

7 7a2/7M75Al 7aiA27a 3 ?0 "A'p'7a4 iV (^)(- T ^ / + 


S3 


= N(P'] 


7 a2 ^Ap7/,75Ai7a 1 A27ay7a 4 7 A A r (^)(V r T - AT) 


5 4 


= N(P r 


7 A '7a 2 «O r A»7/ t 75Al7aiA27a-^0"A'p'7a 4 7 A ^(P)(-TT) 


Si 


= N(P'] 


lajh^il^lcJlc.MP^VV + AA + AV + VA) 


c 


AT ( D'" 

= 1\ 


7 7a 4 ^7^75Al7aiA27a 3 ^A'p'7a2 iV (PJ(-- / 17 + ^ A J 


S3 


= N(P'] 


7 a4 ^A P 7M7 5 Ai7 ai A 2 7ay7 a2 7 A iV(P)(VT - AT) 


S4 


= N(P'] 


7 A '7a 4 ^Ap7M75Al7aiA27a 3 ^A'p'7 a2 7 A ^(^)(-^) 



5! = 7V(P / )7a 2 /7a 1 Ai7M75A27a 3 /7a 4 ^V(P)(W + AA + Ay + yA) 

5 2 = iV(P')7 A '7a^7a 1 Al7/.75A27a 3 ^A'p'7a 4 iV(P) (-T\/ + 774) 

5 3 = iV(P / )7 Q2 ^A P 7 Ql Ai7M75A27ay7 a4 7 A iV(P)(VT-AT) 
<S 4 = iV(P / )7 A '7a 2 ^A P 7a 1 Ai7M75A27a 3 ^v P '7a 4 7 A iV(P)(-TT). 



Appendix C: Results 

We start with the results for G P M and mention that the results for can be obtained 
by the interchange of e u -H- e d . Furthermore, the results for G S A differ from them by the 
absence of the quark charge, the exchange of 7 M — > 7^75, and an opposite overall sign in the 
results of the seventh diagram in both considered quark position cases. 
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First, we consider the case when the discussed quark is at the bottom position. In this 
case, we obtain 

(±^L^L [ M M \(V-A) 2 + 4T 2 ]N(P'h N(P) 



{4nas)2 GU [ , ''I 1 )2 . [d I ] * [(V ~ A? + AT 2 ] N(P> hfl N(P) 

j u 2 {u 2 + u 3 y v 2 {v 2 + v 3 y 



216 g 4 J u 2 {u 2 + u 3 f v 2 (v 2 + v 3 ) 



i^sfe u f [du] [dv] ^^p'^N^ 



216 Q A J u 2 u 3 {u 2 + u 3 ) v 2 v 3 (v! + v 3 ) 
(Aixa s ) 2 e u f [du] [dv 



- , f M , M , \-(V - A) 2 ] N(P'h,N(P) 

(±^L^L j H [ {V + A) 2 ] N(P'h,N(P). 

Second, we consider the case when the discussed quark is at the top position. In this 
case, we obtain 

{A ^ )2ed ' [du] [dv] [2(V 2 + A 2 )] JV(P')7 M iV(P) 



216 Q 4 J Uiim + u 2 ) 2 v^Vi + v 2 f 



(4 ™ s)2ed r [<lll] [<lr] [2{V 2 + A 2 )]N{P') 1 ,N{P) 



216 Q A J u 2 ( Ul + u 2 ) 2 v 2 (v l +v 2 ) 2 



(4 ™ s)2Q I" [du] [dv] [-(V + A) 2 ]N(P'h,N(P) 



216 J uiu 2 (ui + u 2 ) v\v 2 (vi + v 3 ) 
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(4 ™* )2ed f [du] [dv] [-(V-Af]N(P>)^N(P) 



216 Q 4 J uiu 2 (ui + u 2 ) v x v 2 (v 2 + v 3 ) 

(±^l^L [ H IM r 4T 2i N(pf) N(p) 

216 Q 4 7 u 1 u 2 (u 2 + u 3 )v 1 v 2 (v 1 + v 3 ) 1 J 1 J7m 1 J ' 

Finally, we present the results for G\. The discussed quark must be at the bottom 
position here, so that 

(4 ™ s)2 1 f lM [dv] [2(T(V — A) + (V — A)T)\ N(P> h ^N(P) 



216 g 4 7 u 3 (u 2 + u 3 ) 2 v 2 (v 2 + v 3 ) 2 



{AlXas? 1 r h 1 " ] [dv] [2(T(V — A) + (V — A)T)\ N(P'h^N(P) 



216 Q 4 J u 2 (u 2 + u 3 ) 2 v 3 (v 2 + v 3 ) 



(4 ™ s)2 1 f [M [dv] [2(A-V)T]N(P> h ^N(P) 



[2T(A - V)] N(P')^ 5 N(P) 



216 Q 4 J u 2 u 3 {u 2 + u 3 ) v 2 v 3 (vi + v 2 ) 
(4vra s ) 2 1 f [du] [dv] 



216 Q A J u 2 u 3 {u 2 + u 3 ) v 2 v 3 (v! + v 3 ) 
(4tt« s ) 2 1 f [du] [dv 



r [du\ jcb] uy + A)(V + A)]N(P')w 5 N(P). 

216 g 4 7 u 2 u 3 ( Ul + u 2 ) v 2 v 3 ( Vl + v 2 ) Ll A 
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